Dynamic heterogeneities in attractive colloids 
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We study the formation of a colloidal gel by means of Molecular Dynamics simulations of a model 
for colloidal suspensions. A slowing down with gel-like features is observed at low temperatures 
and low volume fractions, due to the formation of persistent structures. We show that at low 
volume fraction the dynamic susceptibility, which describes dynamic heterogeneities, exhibits a large 
plateau, dominated by clusters of long living bonds. At higher volume fraction, where the effect of 
the crowding of the particles starts to be present, it crosses over towards a regime characterized by a 
peak. We introduce a suitable mean cluster size of clusters of monomers connected by "persistent" 
bonds which well describes the dynamic susceptibility. 

PACS numbers; 82.70.Dd, 82.70.Gg, 64.70.Pf 
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Recent advances in colloidal science allow to obtain col- 
loidal particles or nanoparticles with specific functional 
properties, of electronic, chemical, biological or mechan- 
ical nature. Hence the packing and aggregation of col- 
loidal particles are important for a wide variety of appli- 
cations, including biological arrays, sensors, paints, ce- 
ramics, and photonic crystals. 

In particular, by adding polymers to a colloidal so- 
lution, it is possible to induce an effective short range 
attraction between colloidal particles, known as de- 
pletion effect. Attractive colloids exhibit a rich phe- 
nomenology in the temperature - volume fraction plane 
dl, im II IE SB]: At high temperature the attraction 
can be neglected and, by increasing the volume fraction, 
if crystallization is avoided, a hard sphere glass transi- 
tion occurs. By decreasing the temperature, the effect of 
the short range attraction induces an attractive glass line 
which is strongly temperature dependent, as predicted 
by Mode Coupling Theory Q, and by mean field the- 
ory 0- This line at low volume fraction is identified 
as a gelation line. In fact at low temperature and low 
volume fraction, attractive colloids are known to exhibit 
a structural arrest with properties similar to the sol-gel 
transition [1, [l^, [HI, 0, [S] ■ Although intensively stud- 
ied both experimentally and numerically, the theoretical 
understanding of colloidal gelation, compared with chem- 
ical (irreversible) gelation and glass transition, is still far 
from being reached. 

A promising approach to the comprehension of the 
complex dynamics of such systems is based on the study 
of the dynamical heterogeneities. This concept has been 
successfully introduced in glasses [3, EM 16, [l3) IHj 
[igl . |20| . to take into account the correlated motion of 
particle clusters, which decorrelates after a characteris- 
tic time T, of the order of the relaxation time. These 
dynamic heterogeneities are described quantitatively via 



the so-called dynamic susceptibility [l6|. In a Lennard- 
Jones binary mixture, a typical model for glass transi- 
tion ^21|], the dynamic susceptibility grows as a function 
of the time, reaches a maximum and then decreases to 
a constant [20| . consistently with the transient nature 
of the dynamic heterogeneities. Dynamic heterogeneities 
have been also observed in attractive colloidal systems 
in both exp eriments [2^, 113] and numerical simulations 
(25I l26l . lirf . In particular in 27 1 a systematic study of 



the dynamic susceptibility was done along the attrac- 
tive glassy line. Typically the dynamic susceptibility dis- 
plays a well pronounced peak, however in the attraction- 
dominated limit, the dependence on both time and wave 
vector markedly differs from that in standard repulsion- 
dominated systems (hard-sphere limit). 

A rather different behavior was instead found recently 
in a model for irreversible gels, made of monomers with 
permanent bonds [22|. It was shown in fact that the 
dynamic susceptibility, defined as the fiuctuations of the 
self Intermediate Scattering Function (ISF) [l^, in the 
limit of low wave vector, fc ^ 0, tends for long times 
to a plateau, whose value coincides with the mean clus- 
ter size. As a consequence, as the system approaches 
the gel transition (i.e. the percolation threshold), the 
value of the plateau diverges. Such finding shows that, 
in irreversible gelation, the heterogeneities coincide with 
clusters of monomers linked by permanent bonds, and 
differently from glasses, have a static origin and do not 
exhibit any decay. 

In this paper we study the dynamic susceptibility in 
a DLVO type of model for colloidal gelation [28, 29], at 
very low volume fraction and temperature, where the sys- 
tem exhibits more markedly gelation properties. Using 
MD simulations we find an interesting and completely 
unusual behavior. The dynamic susceptibility, for low 
wave vectors, increases with time until it reaches a large 
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plateau value of the order of the mean cluster size, as 
found in the model for irreversible gels, and for time 
larger than the bond lifetime, it decays to 1. Only at 
higher volume fractions, where the bond lifetime is com- 
parable to the density-density structural relaxation time, 
the dynamic susceptibility exhibits a crossover towards 
the glassy behaviour with a well pronounced peak. In- 
terestingly, a proper definition of a time dependent mean 
cluster size, where a cluster is made of "mobile" particles 
connected by bonds present at both time and t, well 
reproduces the observed behavior of the dynamic sus- 
ceptibility for a fixed low temperature and low volume 
fractions. At higher volume fraction, however, the dy- 
namic susceptibility starts instead to display a peak and 
shows a discrepancy with the time dependent mean clus- 
ter size. This indicates a crossover towards a new regime 
where, besides the clusters, also the crowding of the par- 
ticles starts to play a role in the slowing down of the 
dynamics. Although we have considered a DLVO type 
of model, we expect that the main results of our paper 
should be valid also for other model systems exhibiting 
colloidal gelation 

The model - The DLVO model, considered here, has 
been previously studied using Molecular Dynamics (MD) 
In agreement with experimental findings 
this model displays a structural arrest very close 
to the percolation threshold at low temperature, where 
clusters are made of particles connected by long living 
bonds [33|. 

In Refs. [s^l, we found that, at low temperatures, 
increasing the volume fraction, the system undergoes a 
transition from a disordered cluster phase to an ordered 
hexagonal lattice of tubular structures and, at higher vol- 
ume fraction, to an ordered lamellar phase. If this or- 
dered state is avoided, the system enters a "supercooled" 
metastable liquid phase until structural arrest (gel) oc- 
curs 13, 13] very close to the percolation threshold. For 
this reason we introduce here a small degree of polydis- 
persity, which actually hinders the formation of the or- 
dered phases. In this way we are able to fully investigate 
and characterize the gel formation at different temper- 
atures and volume fractions. The interaction potential 
between two particles i and j is: 




(1) 



75.08, C ^ 0.49 



where A = 3.56, B = 7.67, C = 
[35| . and aij = (cr^ + o"j)/2. In the following we con- 
sider the radii ct^ randomly distributed in the interval 
a - S/2 < ai < a + S/2 with S = 0.05cr. The potential 
is truncated and shifted to zero at a distance of 3.5(7. 
The temperature T is in units of e/ks, where ks is the 
Boltzmann constant, wavevectors are in units of cr"^, and 
times in units of ma'^/e, where m is the mass of the 
particles. The volume of the simulation box is kept con- 
stant, V = 50007rcr'^/3, and different volume fractions are 
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FIG. 1: (color online). The self ISF, Fs{k,t), for T = 0.15, 
fcmm = 2-!v/L ~ 0.36, and (p = 0.01, 0.05, 0.08, 0.10, 0.11, 
0.12, 0.13 (from left to right). 



obtained by varying the number of particles N, so that 
(f) — na^N/GL^. We have performed Newtonian MD at 
constant NVT using the velocity Verlet algorithm and 
the Nose- Hoover thermostat with time step At — 0.01. 

At low temperatures and low volume fractions, we find 
the same properties as in the monodisperse case [l^l: A 
cluster phase followed by a percolation transition in the 
same universality class as random percolation. 

Self Intermediate Scattering Function- In order to 
study the dynamic behaviour of the model, we measure 
the self ISF, F,(fc,i) = ^($^(fc,i)), where $,(fc,i) = 



N „jfe-(n(t)-n(o)) 



. .) is the thermal average, and 



ri{t) is the position of the i-th particle at the time t. 
In Fig.[T]Fs(fc, t) is plotted for T = 0.15 and wave vector 
kmin — 'i.TT / L ~ 0.36. Increasing the volume fraction, 
the relaxation functions show a two-step decay, with a 
plateau value decreasing with the wave vector; The long 
time tail is fitted by stretched exponentials. 

From the self ISF we calculate the structural relax- 
ation time, Ta{k), defined as Fs{k,Ta{k)) ~ 0.1. In 
Fig. [2] Ta{kmin) IS plotted as a function of the volume 
fraction for two values of temperature, T = 0.15 and 
0.25. In the same figure for the sake of comparison 
we also plot the bond lifetime, Tf, [s^, defined through 
the bond correlation function as B{ti,) ~ 0.1, where 

s:i,[("M(t)"ij(o))-<"ij)' ' 



B{t) 



-, nij{t) = 1 if particles i 



and j are bonded at time t, (t) = otherwise [33| . The 
data show that, at the lower temperature, T = 0.15, and 
at low volume fractions, where the bond lifetime is larger 
than the structural relaxation time, Ta{kmin) can be fit- 
ted by a power law, (^c — cj))~°'. In this regime the system 
dynamically behaves as if made of permanent clusters as 
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FIG. 2: (color online). The structural relaxation time, 
Taikmin) (circles) , Compared with the bond relaxation time, 
Tf, (stars), for T = 0.15 and 0.25 (from bottom to top). The 
continuous line is a power law fit (0.14 — (j!>)~^'*. 



in irreversible gels. At the higher temperature, T = 0.25, 
instead the bond lifetime is of the same order of magni- 
tude than the structural relaxation time, and we observe 
a smooth increase of TQ(fc„jj„) not associated to a power 
law critical behavior. 

Dynamic heterogeneities - According to the above in- 
terpretation of the dynamics, we should expect, for time 
windows comparable or larger than Tf,, a deviation from 
the dynamics of irreversible gels even for T — 0.15. This 
crossover in the dynamics can be better stressed by look- 
ing at the dynamic susceptibility, defined as X4(fci^) — 

i[(|$,(A;,i)|2}-|(ci>,(fc,i))p]. 

In the main frame of Fig. [31 x^^i^mm^t), is plotted as 
function of t, for T — 0.15 and different (f). For low 
volume fractions, after a time roughly of the order Tq, 
Xiif^mimt) reaches a plateau, close to the value of the 
mean cluster size, S = s^Us/ ''^s (see also the inset 
of Fig. [3] where the maximum of the dynamic susceptibil- 
ity is compared with the mean cluster size), and decays 
after a time roughly of the order t\,. Only for times much 
smaller than rt, the system behaves as if the bonds were 
permanent. As the volume fraction increases, and Tq, 
become of the same order of magnitude, the plateau dis- 
appears and Xii^min, t) exhibits a sharp peak similar to 
those of glassy systems. We have also checked that for 
temperature T = 0.25, Xii^mimt) exhibits a small peak 
and no plateau. 

The long time behaviour of Xi{km.i7i,t) can be quan- 
titatively described in terms of a time dependent mean 
cluster size, S{t). The function S{t) is the mean cluster 
size defined by persistent bonds, that is bonds that are 



FIG. 3: (color online). Main frame: The fluctuations of the 
self ISF, X4(fc™in,t), for T = 0.15 and = 0.01, 0.05, 0.08, 
0.10, 0.11, 0.12 (from left to right). Inset: The mean cluster 
size, S (triangles), as a function of the volume fraction, 0, 
compared with the maximum of the dynamic susceptibility, 
X4,{kmin,t*) (circles). 

present at both time and time t [s^] (see inset of Fig. d]). 
To describe also the short time behaviour of Xii^^mimt) 
we note that the main contribution to Xiif^mimt) comes 
from the "mobile" particles [38]. We therefore modify 
further the definition of the mean cluster size introducing 
a mean cluster size for "mobile" particles, S„i{t), where 
the clusters are made as before by particles connected 
by persistent bonds, but restricted to particles that in 
the time interval [0, t] have moved at least a distance tq. 
For each volume fraction, we fix ro so that the time at 
which Sm(t) and Xiikmimt) start to grow is the same. 
We find that tq depends on the volume fraction roughly 
as In main frame of Fig. IH Smit) is compared 

with Xiif^min , t) . At least at very low volume fractions, 
the two quantities agree surprisingly well. At = 0.10 
a deviation begins to appear, that becomes manifest at 
(j> = 0.12, in the shape of a peak that grows to values 
higher than Sm(t) or S{t). This peak, which appears 
when Tq and Tf, are of the same order of magnitude shows 
that the contribution to dynamic correlations, comes not 
only from the presence of long living clusters but also 
from the crowding of the particles. 

In conclusions we have shown that in a model for col- 
loidal systems at low volume fraction when the bond life- 
time is much smaller than the relaxation time, the dy- 
namic heterogeneities can be described in terms of time 
dependent clusters made of "persistent" bonds. We sug- 
gest that the analysis presented here should apply not 
only to colloidal gelation at low volume fraction, but 
also to micellar system, where a crossover from gel-like 
to glass-like behavior was found experimentally and nu- 
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FIG. 4: (color online). Main frame: The dynamic suscepti- 
bility, X4(fcmm,i) (circles), compared with the time depend- 
ing mean cluster size of mobile particles, Sm{t) (stars) for 
r= 0.15 and 0.05, 0.10, 0.12. Inset: x4kmin,t) (cir- 
cles) compared with the time depending mean cluster size, 
S(t) (stars) for the same temperature and volume fractions. 

merically [41]. These findings offer a unique and coherent 
interpretation of dynamic heterogeneities in gelhng sys- 
tems. Finally, we also suggest that the concept of time 
dependent cluster, considered here, may be generalized to 
give a satisfactory definition of dynamic heterogeneities 
also for Lennard- Jones or hard sphere glasses. 
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